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An effective field theory is used to describe light nuclei, calculated from quantum chromodynamics 
on a lattice at unphysically large pion masses. The theory is calibrated at leading order to two 
available data sets on two- and three-body nuclei for two pion masses. At those pion masses 
we predict the quartet and doublet neutron-deuteron scattering lengths, and the alpha-particle 
binding energy. For m,r=510 MeV we obtain, respectively, "^OnD = 2.3 ±1.3 fm, ^Ono = 2.2 ±2.1 fm, 
and = 35 ± 22 MeV, while for m.n-=805 MeV "^ano = 1.6 ± 1.3 fm, ^OnD = 0.62 ± 1.0 fm, and 
Ba = 94 ± 45 MeV are found. Phillips- and Tjon-like correlations to the triton binding energy are 
established. Higher-order effects on the respective correlation bands are found insensitive to the pion 
mass. As a benchmark, we present results for the physical pion mass, using experimental two-body 
scattering lengths and the triton binding energy as input. Hints of subtle changes in the structure 
of the triton and alpha particle are discussed. 
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I. INTRODUCTION 

The vast number of phenomena of the nuclear chart depend on a relatively small set of quantum chromodynamics 
(QCD) parameters — in the low energies relevant for nuclear physics, a mass scale Mqcd associated to the strong 
coupling constant, the masses niq of the two lightest quarks, the electromagnetic coupling strength, and the vacuum 
angle. Lattice QCD (LQCD) is a numerical framework which enables us, at least in principle, to relate nuclear and 
QCD parameters, once effects due to finite lattice spacing and size are removed. The last few years have witnessed 
significant progress in predicting the properties of light nuclei with nucleon number A < 4, but at relatively large 
quark masses and neglecting time-reversal and isospin violation. (See Ref. [T] for a review and a list of relevant 
references). 

Increasing A at fixed quark masses presents significant difficulties because the noise-to-signal rate increases expo¬ 
nentially. Although there seem to be ways around this problem [T], large A also requires that longer distances be 
covered by the lattice, since the nuclear volume increases with A. As in other areas of physics, it is profitable to 
change to a more effective description, in this case to an effective held theory (EFT) involving nucleons as degrees of 
freedom. Because an EFT is based on the most general Hamiltonian with the appropriate symmetries, it is guaranteed 
to produce ^-matrix elements consistent with the S matrix of the underlying theory [2, here QCD. After matching 
the EFT amplitudes to the LQCD-calculated quantities at small A, one can describe the longer-distance dynamics 
involved in larger-A systems within the EFT [5], which is considerably simpler than doing so directly within LQCD. 

Most LQCD results so far concern binding energies, but reactions convey much more information and will command 
increasing attention in the years to come. Unfortunately, as discussed in Ref. [3] , which also summarizes the progress 
in this held, volume artifacts are more pronounced. EFT quite naturally accounts for scattering states, and allows 
bound states and scattering to be treated on equal footing. Here we elaborate on the hndings of Ref. [2] for A < 4 
and extend, for the hrst time, LQCD predictions to reactions involving nuclei. As an example, we consider neutron- 
deuteron (nD) scattering at low energies, where the two 5'-wave channels — with total spin s = 3/2 (quartet) and 
s = 1/2 (doublet) — are most important. 

The noise-to-signal rate in LQCD also increases with decreasing mq. Results obtained with unphysical rriq can, in 
principle, be extrapolated to the physical point in a systematic way using chiral effective held theory (y^EFT), as long 
as pion masses are within the radius of convergence of the latter. From yEFT with up to one nucleon — that is, 
chiral perturbation theory (xPT) — one obtains the niq dependence of, for example, the average pion mass (wTr) [H], 
and of the nucleon (m^) and Delta (ttia) masses [6]. The niq dependence of some few-nucleon observables has also 
been estimated |3, but unfortunately signihcant uncertainties still exist due to subtleties in the proper accounting of 
renormalization-group (RG) invariance in this non-perturbative context [71 1561 IS5] . 

The average pion mass is commonly used as a measure for the detuned value of the average quark mass. At 
present, LQCD can be carried out in the meson and single-hadron sector down to values of close to physical, 
where the low-lying mass spectrum is reproduced within theoretical error bars (see Ref. |12j for a status report). 
Comparison with LQCD data suggests that xPT converges for pion masses no larger than about 500 MeV [53]. In 
contrast, the quark masses employed in current nuclear LQCD are likely beyond reach of xEFT. 

As proposed in Ref. [5] and elaborated here, the EFT that describes existing light-nuclear LQCD data need not 
include pions explicitly. In fact, it has been understood for over fifteen years that even at the physical pion mass light 
nuclei are well described by pionless EFT (•/EFT), an EFT with non-relativistic nucleons interacting through contact 
forces with an increasing number of derivatives — each with a strength parameter or “low-energy constant” (LEC) 
— which contribute at increasing orders. In two-nucleon scattering, ■/'EFT reproduces [131 f55H37] the effective range 
expansion (ERE): scattering lengths at leading order (LO), effective ranges at next-to-leading order (NLO), etc. It 
thus also gives two-nucleon binding momenta in the and ^Sq with corresponding accuracy. More importantly, 
•/EFT offers a consistent extension of the ERE to other systems [30] . For example, S'-wave nD scattering in the quartet 
channel can be very accurately postdicted [101 - 13(1] once the two-nucleon LECs have been fixed in the two-nucleon 
system. In the doublet channel, in contrast, RG invariance requires that the three-body force with no derivatives 
appear already at LO, with isospin-symmetric corrections starting beyond NLO [MllIIHMj. Gurrent evidence from 
the RG in the four-body system suggests that there is no four-body force up to NLO [30M3] . The existence of a single 
three-body parameter up to NLO, which determines the three- and higher-body spectra, leads to many correlations 
among few-body observables at fixed two-body input. Examples are the so-called Phillips [33] and Tjon m lines 
obtained in plots of the doublet nD scattering length [331131] and alpha-particle binding energy [331131] as functions 
of the triton binding energy. Higher partial waves in three-nucleon scattering [11131, four-nucleon scattering [36] . 
and even ®Li [33] can also be reasonably well described in /-EFT. 

We will show that an analogous approach to describe light nuclei is equally useful at larger mT^. Using higher-than- 
physical niq not only increases TOtt, but also changes the nucleon mass ttin and the masses of all other hadrons. We 
will argue on the basis of scales inferred from LQGD data that nucleons are sufficient for momenta up to with 
neither explicit pions nor other baryons. Whether it is indeed (instead of, say, toa — Wn) that determines the 
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convergence rate of the theory used here is the subject of an upcoming investigation. At each value of rriT^ a pionless 
EFT exists with specific values of the LECs; we refer to ^fEET with varying as X'EFT in the following. Until 
nuclear LQCD calculations are extended to include time-reversal and isospin violation, is the only QCD parameter 
determining nuclear properties. Existing data at = 805 m and 510 [55] MeV give A < 4 binding energies that 
are much larger than in the real world and increase with the pion mass. The dineutron is bound, which could signal 
qualitative new features in lattice worlds. An obvious question is the extent to which properties of ^fEET survive in 
yf'EFT, where all scales change. 

In Ref. |2] the binding energies of nuclei with A < 6 were studied in LO yf'EFT using as input the LQCD data 
for dineutron, deuteron and triton/helion at = 805 MeV |37j . The alpha-particle binding energy provided a 
consistency check between X'EFT and LQCD data, and the A = 5, 6 binding energies obtained with X'EFT can be 
viewed as an extrapolation of LQCD. Here, we extend X'EFT to the = 510 MeV LQCD data |55| and to a broader 
range of observables including scattering amplitudes. 

The methods of ^EFT have for some time been deployed in the study of reactions directly on the lattice [H HD] . 
Both two-nucleon elastic scattering |39j and neutron radiative capture on the proton im have been considered directly 
on the lattice. Our strategy is, instead, to analyze reactions outside the lattice box with X'EFT once its LECs have 
been determined from binding energies at LO and, eventually, also two-nucleon scattering observables at NLO. We 
exploit the dramatic advances in the development of the so-called ab initio methods that have taken place in nuclear 
physics over the same period in which ^EFT was formulated. In particular, here we employ the effective-interaction 
hyperspherical-harmonic (EIHH) method |42ll44j . and the refined resonating-group (MGM) method |69|. Although 
these methods have been developed for traditional nuclear potentials, they can be adapted to pionless EFT, as already 
done for X'EFT in the EIHH |5] and /EFT in the KGM [55] . 

Thus, we show that /EFT remains useful in nuclear systems with A < 4 and extrapolate LQGD data to observables 
that might not be as easily obtained in the lattice. This is analogous to the use of /EFT correlations [13 SH] to infer 
values of poorly measured observables in the real world. If and when scattering observables are determined directly 
on the lattice, our predictions will be a further test of the consistency between XEFT and LQGD, establishing the 
validity of a theory with only contact interactions over a range of TOtt from 140 MeV up to 805 MeV. Such a consistency 
would provide a benchmark for the extension of this method to the less-understood xEFT, once LQCD data reaches 
sufficiently small pion masses. 

We summarize the article as follows. In Sec. jTTjwe discuss the degrees of freedom and breakdown scale ofXEFT for 
up to ~ 800 MeV. Still in Sec. jH] we present the LO Hamiltonian and the regulator we use in calculations, which 
employ the computational tools introduced in Sec. 


HI the EIHH and MGM methods. In Sec. IV 


we determine 

the LO LECs from the LQCD data for A < 3 in the alternate reality assessed via LQCD at various m.,^. With 
the Hamiltonian thus calibrated, we calculate in Sec. jV] the alpha-particle binding energy, establish the heavy pion 
Phillips and Tjon lin es, a nd predict the doublet and quartet neutron-deuteron scattering lengths ^UnD and “^anD- As 
we conclude in Sec. VI the procedure is analogous to the development of /EFT over the last decades, namely, a 


calibration of a small set of parameters to data in order, first, to obtain predictions of low-energy observables and, 
second, to explain correlations amongst them. 


II. PIONLESS EFFECTIVE FIELD THEORY 

At physical 77i,r, pionless EFT with nucleons as the sole degrees of freedom has proved useful for light nuclei in the 
low-momentum regime — see Refs. |43|48] for reviews and Ref. [49] for a pedagogical introduction. Its organizational 
scheme (“power counting”) is based on two basic scales: the breakdown scale Qhigh estimated as and an unnaturally 
small scale H related to the inverse of the two-nucleon scattering lengths ^’^onn in the singlet/triplet channels. 

For external relative momenta k < m.,^/2, the deuteron and the virtual singlet bound state are the only singularities 
of the two-body scattering amplitude. All mesons and excited baryons are short-range effects. The amplitude from 
a Lagrangian built of derivative contact operators made of nucleon fields can be matched to all orders of the ERE. 
Matching the LO amplitude to the ERE results in a Cs,t oc 47r(^’^aNN)/niN ~ 47r/(mNH) scaling of the non-derivative 
four-nucleon contact term. As the scattering lengths ^’^onn are large relative to the pion range l/m^r for both NN 
S'-wave channels, a refined power counting is required |65L 166] that goes beyond naive dimensional analysis. Of 
course, care has to be taken that the necessary regularization and the inclusion of higher-order contributions do not 
introduce poles within the radius of convergence. As long as those singularities are beyond the pion threshold, m^, 
/EFT converges for two-nucleon processes at momenta Q < Qhighj including the ^Si (deuteron) and poles [T51I57] . 
In LO the two LECs Cs,t suffice. 

Extending /EFT to systems with more nucleons requires understanding how H enters the LECs of multinucleon 
interactions. The fact that the non-derivative six-nucleon contact interaction is needed to define the EFT at LO 
[mug implies its LEG scales as Dd (47r)^/(mNH^). In contrast, the apparent lack of similar RG enhancements in 
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TABLE I. Relevant scales for a low-energy nuclear effective field theory. Physical data in the first column is relevant for /fEFT, 
lattice data summarized in the second and third for/PEFT. All numbers are given in MeV. 


pion mass 

771 TT 

139.5 ±0.1 [72] 

511 ±2 [38] 

806 ± 1 [37] 

nucleon mass 

mN 

939 ± 1,938 ± 1 [Z3] 

1320 ± 3 [38] 

1634 ± 18 [37] 

Delta-nucleon mass difference 

(5a = niA — mN 

292 ± 1 [50] 

« 200 [SI] 

« 180 [SI] 

dineutron binding energy 


— 

7.4±2 [3S] 

15.9 ±4 [37] 

deuteron binding energy 

Bn 

2.22 [52] 

11.5 ±2 [38] 

19.5 ±5 [37] 

triton binding energy 


8.482 [73] 

20.3 ±4.5 [38] 

53.9 ± 10.7 [32] 

inverse singlet scattering length 

^np 

-8.31 [53] 

n.a. 

84.7 ± 18 [39] 

inverse triplet scattering length 

^np 

36.4 [53] 

n.a. 

108 ± 13 [H] 

Delta effective momentum 

'\/2771i\r^A 

741 

890 

767 

two-nucleon binding momentum 

TTIn (-Snn + ) /2 

46 

112 

170 

triton-to-deuteron binding ratio 

St/-Bd 

3.82 

1.8 

2.8 


other contact interactions suggests they appear only in higher orders. 

As -^EFT is applied beyond the deuteron, one needs to account for effects of the Coulomb force among protons. The 
importance of Coulomb effects is characterized by a ratio arriN/Q, where a is the fine-structure constant. Although 
crucial for very low-energy proton-proton m and proton-nucleus scattering, the Coulomb interaction should be 
subleading in relatively deep ground states such as helion and alpha particle, where binding momenta are much larger 
than aruM- 

At LO, the /EFT Lagrangian can be written as 

C^O = (ido + — ) N + — (Af^crVV^ N)^ N) + — (N'^N) 

\ 2?7t.n is 8 

+Dd{N^N){N^N){N^N), ( 1 ) 

where is a bi-spinor in both spin and isospin spaces, and cr* (t“) are the spin (isospin) Pauli matrices, the index 
i (a) running over spin (isospin) vector components of the projection operators on the spin singlet (triplet) state. 
Higher orders contain terms with more derivatives and/or nucleon fields, including those necessary to ensure Lorentz 
invariance (in a Q/m^ expansion). 

Somewhat surprisingly, /EFT seems to converge for triton and helion [231 HH US] , and even for the more-bound 
alpha particle |29l|32]. At the physical point, /EFT is useful even at LO to explain features like correlations amongst 
three-body observables (the Phillips line) and between three- and four-body data (the Tjon line), with just the 
neutron-proton scattering lengths as input. With an additional condition which conventionally fixes either 

the triton binding energy or the neutron-deuteron doublet scattering length ^OnD, a few four-nucleon observables 
— e.g., the binding energy of ^He [291132], and the neutron-triton and proton-helion scattering lengths [36] — have 
been found to agree with data within the expected uncertainty margin. The only exception so far seems to be the 
resonance location in the 0“ neutron-triton channel, which was found to be cutoff, and thus renormalization-scheme 
dependent [ST]. The origin of this pathology is unknown. LO results for ®Li [SS] do not allow conclusions about the 
range in A where /EFT converges. 

With the usefulness at physical thus established, we follow an analogous approach at heavier niq. Available 
lattice data m identifies the pion, still, as the lightest meson and the Delta as the lowest excited state of the nucleon. 
However, the ratios between the nucleon, pion, and Delta masses change, see Table |l] Also, nuclei become increasingly 
more bound. 

The relevant momentum is very clear in the two-nucleon system, from either the inverse scattering lengths or the 
two-nucleon binding momentum estimated from the average two-nucleon binding energy. At all pion masses it is much 
smaller than the nucleon mass, meaning nucleons are nonrelativistic, and even than the pion mass itself, ensuring 
pions can be integrated out. However, in contrast to the physical world, m,r > "m-A — Wn = Sa for the two lattice 
simulations, and hence one might wonder if the Delta should not be included as an explicit degree of freedom. 

The reason the Delta can still be integrated out is, of course, that in a nonrelativistic theory the relevant quantity 
for convergence is momentum, not heavy-particle mass. In this case, it is the “Delta effective momentum” •\/2mN^Aj 
which remains above, or at least near, the pion mass. That y/2m^SA is the relevant scale was shown explicitly in 
Ref. [53] for the two-nucleon ^5'o channel. In this case, the lowest accessible state with excitations has two Deltas, 
and in addition to Cg two other non-derivative contact interactions need to be included: two-nucleon/two-Delta and 
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four-Delta. Under the assumption that all three LECs are of a similar size Cq, they scale as [51] 


Cn = 


Itt 


('«NN + '«(234 MeV)-" 


( 2 ) 


at 771^=140 MeV. Because the inverse value of the singlet scattering length for m^=805 MeV, displayed in Table [T] 
is about 10 times larger in magnitude than its physical analog, the ensuing size of Cq would decrease and allow for 
higher typical momenta in the two-nucleon amplitude. However, the Delta effective momentum is still several times 
larger than the inverse scattering length. Operators in a Deltaful, pionless theory should then show similar scaling 
behavior as for physical pion mass, where the Delta can be integrated out. Removing the Delta generates an effective 
range not accounted for in LO, but this contribution is characterized by the Delta effective momentum, which does 
not seem to be smaller than the inverse pion mass. 

This argument can be generalized to other channels where ^/2rrh^SA is replaced by -x/rr^^A^, with the 
difference between the mass of the state containing the nucleon excitation(s) and 2mN. Since lattice results suggest 
that the lowest state with a single excitation involves the Roper resonance, whose mass is somewhat larger than the 
Delta, one does not expect a significant decrease in convergence rate by keeping only the nucleon explicit in the EFT. 

Therefore, we formulate yPEFT as an EFT formally equivalent to -jfEFT, but with different scales and values for 
the LECs. The breakdown or high-momentum scale Qhigh is assumed to be the smallest of and The 

low-momentum scale Qiow is set by the binding momenta of the nuclei we consider and by the external momenta in 
the reactions we are interested in. We expand all observable in powers of Qiow/Qhigh- Eventually, an NLO calculation 
will yield an estimate on the convergence rate and thereby the breakdown scale of yf'EFT. The Lagrangian in LO is 
given by Eq. Q, in which four m 7 r-dependent parameters enter: the nucleon mass Ton and the LECs Cg^t and Dd- 

For the calculation of few-body observables we solve the Schrodinger equation in configuration space. The potential 
is the sum of all irreducible contributions to the H-body scattering matrix from the Lagrangian. This amounts 
at LO to the sum of three tree-level diagrams with vertex factors Cg^ and Dd- The infinities resulting from the 


zero-range contact interactions are here regularized via Gaussian regulator functions, exp(—A^r?j/4)/(167r^/^) for 
two nucleons i,j and A® exp[—A^(r^^ -|- r^j;,)/4]/(647r^) for three nucleons i,j, fc, where A arbitrarily separates states 
included explicitly as propagating degrees of freedom from states accounted for implicitly in the LECs. If it is smaller 
than the breakdown scale it produces larger errors than the truncation of the EFT Lagrangian (Sec. jV A^ exemplifies 
ramification of a violation of this condition). The resulting Schrodinger equation for the A-body wave function and 
the corresponding energy E takes the form 


— —A— {Cl— Cq^i) (Ti ■ (Tj\e ’’o-I- i Cik+^jk) K (3) 


2m^ 

i<j i<j<k eye 

Here, a factor from the regulator was absorbed into the bare couplings of Eq. Q: 

A3 


^ 04 (A) = ^^C.(A), 


A3 


^Lo(A) = Y^a(A), 


(4) 

(5) 

( 6 ) 


As in any EFT, the bare LECs depend on A so as to guarantee that observables do not. The A-dependent LECs are 
determined from input data in Sec. IV after we discuss the solution of Eq. ([^ in the next section. 


III. TOOLBOX 

To solve the Schrodinger equation we have utilized two computational methods: EIHH and KGM. Hereafter, we 
present a short description of both methods. 


A. The Effective-Interaction Hyperspherical Harmonics Method 

The hyperspherical coordinates are the D-dimensional generalization of the 3-dimensional spherical or polar co¬ 
ordinates. As such they allow the description of the A-body wave function in terms of a single length variable, the 
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/lyper-radius p, and {D — 1) /lyper-angular variables n |75l I76j . Removing the center-of-mass coordinate, the A-body 
dynamics can be described by 7l — 1 Jacobi vectors rji, therefore D = 3A — 3. 

The nice feature of these coordinates is that, in perfect analogy to the two-particle case, the kinetic energy operator T 
of the A-particle system splits into a hj/per-radial and /lyper-centrifugal terms, with a hyper spherical angular momentum 
operator K that depends on Q. The resulting A-particle Hamiltonian reads 


where Ap is the hyper-radial Laplacian. 

The hyperspherical harmonics (HH) y[K] are the H-body generalization of the spherical harmonics. As such they are 
the eigenfunctions of with eigenvalues K{K -I-3A — 5). They form a complete set of hyper-angular basis functions. 
Choosing a complementary set of hyper-radial basis states i?„(p), the A-body wave function can be expanded in the 
form 


^ C^[K]Rnip)yiK]m ( 8 ) 

nlK] 

with coefficients Cn[K]- The nuclear wave function must be complemented by the spin-isospin parts, and the whole 
function must be antisymmetric. The construction of antisymmetric HH spin-isospin basis states is a non-trivial task, 
which, however, has been solved in Refs. [57l[58] . 

To accelerate the convergence rate of the HH expansion, Eq. Q, we construct an effective interaction (El) using 
the Lee-Suzuki similarity transformation m Applying the this method to the HH basis we identify the model space 
P with all the A-body HH states such that K < Kmax, and the complementary space Q = 1 — P as the rest of 
the Hilbert space. The Lee-Suzuki method then gives a recipe to construct a similarity transformation such that the 
spectrum of the resulting effective P-space Hamiltonian, iJ = T -\- coincides with the spectrum of 

Finding however, is as difficult as solving the original problem, and therefore we do not search for the total 

El, but for a partial El constructed through the solution of the simpler two- and three-body problems. 

The resulting El is tailored to our HH model space, and constrained to coincide with the bare one when enlarging 
the model space. This EIHH method has been successfully applied to the study of bound states and reactions 

for nuclear systems with 3 < A < 7. 


B. The Refined Resonating-Group Method 

In contrast to the EIHH method where the few-body wave function is expanded over a complete set of states, the 
KGM is a variational approach that utilizes an over-complete set of states (for its original formulation, see Ref. [7IJll71j : 
for the refinement and implementation, Ref. [53]). To construct these states, the KGM considers all possible channels 
{[c]}, where each channel consists of a specific spin-isospin configuration S[c], a set of Jacobi vectors rji,... 
and the angular momentum quantum numbers fi,... ,£a-i associated with these vectors. The orbital functions are 
given by the ansatz 


A-l 

Rule] iVi,---, Va-i) = n , (9) 

i=i 

where Yim are the spherical harmonics, and are a set of width parameters used to expand the wave function, i.e., 
the sum over the channels includes an expansion of each radial dependence in Gaussians with widths {unj}- 

The few-body wave function is then a linear combination of an antisymmetric product of a spin-isospin channel 
state and the orbital function, coupled to yield the desired total angular momentum quantum numbers JM, 

'I'jM = A'^Cnlc] [i?„[c]G)S[c]]'^^^ . (10) 

n[c] 


The sum over channels allows the consideration of all possible spin-isospin configurations or clusters Sj^]. In practice, 
however, our implementation omits channels that have negligible contribution to the wave function. For example, the 
ansatz for the alpha-particle wave function includes triton-proton, helion-neutron, and deuteron-deuteron spin-isospin 
configurations. The conceivable two-neutron-two-proton arrangement was found to contribute less than 100 keV to 
Hq, and therefore is not included in the variational ansatz. 
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Thus, the KGM method includes three intertwined expansions: i) the cluster or resonating-group expansion, 
which defines the spin-isospin conhguration and the Jacobi coordinates; ii) the partial-wave expansion; and Hi) the 
expansion in Gaussian functions. Convergence is assessed along each of those “axes”. First, the thresholds of a system 
serve as guidance for the initial choice of resonating groups. Second, contributions from subleading partial waves are 
considered. For s-wave nuclei, and central forces, i > 0 configurations do not have to be included due to the cluster 
expansion. Consequently, at this order of our EFT we consider only £ = 0 terms in our description of the light, A < 4 
nuclei. Third, the set of Gaussians is extended and scaled until this modification of the model space does not affect 
binding energies by more than 1%. 

With the KGM, we also calculate scattering observables. To solve the few-body problem with the KGM for a range 
of cutoff (A as introduced above to obtain the regularized Eq. values, i.e., to approximate a wave function with 
structure around rij « A“^, the variational basis has to be either very large — leading to numerical instabilities— or 
tailored to each A — requiring a convergence check with regards to all parameters of the basis set. 

Our variational approach is analog to Kohn-Hulthen’s method m which minimizes a functional parameterizing 
the reactance matrix, corresponding to Ricatti-Bessel asymptotic solutions for uncharged particles and Coulomb 
functions for charged fragments. We use in- and out-going waves as boundary conditions (spherical Hankel functions 
h^), because this method turned out to be more accurate in practice. Specifically, for two-fragment scattering with 
an incoming channel c we denote the relative intercluster Jacobi coordinate by r]c and make the ansatz 


4/ =Ai-h^{T]c)+'^Scc'h+,{T]c>)+ ^ 


[c] 


( 11 ) 


n[c]GC 


with variational parameters Sec' (the S matrix) and If either target or projectile are compound objects, e.g., 

the deuteron in Sec. V A their wave functions are predetermined via the ansatz in Eq. (10) and multiplied with the 
asymptotic solutions of the relative motion. For small distances t]c, the interaction between nucleons of different 
fragments is non-zero and the full scattering wave function will differ from the asymptotic form as given by the first 
two terms in Eq. (111. This difference is described by the third term in Eq. 0- Convergence and stability are 


assessed with respect to the subset C which is taken from the full set of channels. It is sufficient to include those n[c] 
in C which are non-zero for rje ~ A“^ and, as Gaussians, square-integrable. For rjc ^ A“^, this expansion should be 
zero, ie., is identical to the asymptotic solution. The Kohn-Hulthen variational condition expressed in terms of 
the scattering matrix is 




( 12 ) 


where „ is the center-of-momentum energy. This condition yields optimal values for Sec' and iJ„[c]. Here the 
channel index c discriminates between different two-body fragmentations, e.g., neutron/deuteron or neutron/neutron- 
proton singlet, and angular momentum. Using an appropriate decomposition of the Hamiltonian (for the latest 
summary and references to the original work see Ref. [H2), the variational coefficients Sec'^ Dn[c] can be expressed 
in terms of integrals of the short-ranged part of the potential. Therefore, an accurate expansion of the asymptotic 
solution is required for a finite range. In practice, we minimize 





iv) / . l^rafc] l^wfc] 

n[c] 


rje 


(13) 


to approximate the Hankel or Coulomb functions. Finally, we obtain scattering lengths from the phase shift ) 

at a finite through 


a(E_,) 


1 

kcot5{E^„,) 


(14) 


As the scattering length is defined for E^ ^ = 0, the uncertainty due to this approximation has to be assessed. In 
this work, we extracted a at 0.001 MeV, used 10-13 Gaussians to expand the deuteron and singlet neutron-proton 
fragment in the three-body scattering calculations, and adapted the Hankel functions with a weight e = 0.03 fm~^. 
To conclude, we summarize the convergence check: 


• First, we determine appropriate Gaussian basis for the fragments by hxing the number of Gaussians and optimize 
their widths via a genetic algorithm . 









A’moi Kmaa: 

FIG. 1. Dependence of the alpha-particle binding energy (in MeV) calculated with the EIHH method on the maximal 
/ij/per-angular eigenvalue Kmax- Results are shown for a^EFT interaction with A = 2 fm“^ at m.„= 140 MeV (left panel), 
510 MeV (center), and 805 MeV (right). The horizontal green line represents the corresponding RGM result. 


Second, we diagonalize the Hamiltonian, Eq. [Sjin the scattering basis. This basis uses a different coupling 
scheme which adopts the one implied in Eq. (10) for each fragment. The total fragment spins are coupled 
to a channel spin which forms, with the orbital angular momentum on the relative coordinate r]c between the 
fragments, the total J. We enlarge C until the lowest eigenvalues reproduce the thresholds defined by the ground 
states of the fragments and the bound states of the compound system of the two fragments, if there is a bound 
state in the channel (the triton in Sec. VA). 


Third, we take 


lim 1(e) and lim a{E) (15) 

e-fO E—^0 


in Eqs. (13 and ( [T4| . While taking both limits, we identify a plateau in the predicted scattering length for 
e < A2 and < 0.0001 MeV. 


After these steps, we deem the basis large enough for an accuracy that is then dominated by the higher-order 
contributions of the EET expansion. 


C. Comparison 

With EET parameters calibrated as described below, we compared the results for and of the KGM with the 
corresponding EIHH values to test the accuracy of the resonating-group expansion. As an example, we show in Fig. 
the convergence of EIHH calculations to the KGM results for B^ at a cutoff A = 2 fm“^. For all three pion masses, 
the EIHH converges with Kmax to the respective KGM value. 

For subsequent calculations the KGM was chosen to minimize computing time. 


IV. CALIBRATION 

Through Eq. all LO predictions depend on three LECs Cs,t {S,T = 0,1 or 1,0) and Di, besides the nucleon 
mass. Lattice data are available for the two- and three-body binding energies {B^^, B^, B.^.) at TO,r= 510 MeV [38] 
and at m,r= 805 MeV m- At 805 MeV [31], the singlet and triplet scattering lengths (^Onp,^ Unp) and effective 
ranges are also available . We fit Cs,t to the two-nucleon binding energies in the singlet and triplet channels (i?^„, B-^) 
by solving the Schrbdinger equation via the Numerov algorithm. The Di term is fixed through B^ using the KGM. For 
comparison, we also consider the physical pion mass, m7r=140 MeV, where we fit the experimental singlet scattering 
length, in addition to experimental deuteron and triton binding energies. 

These renormalization conditions determine the A dependence of the LECs. The values of the bare LECs for cutoffs 
A = 2, 4, 6, 8 fm“^ are given in Table [nj and a graphical representation of the fit results is given in Fig. The input 
to calibrate the values at the physical pion mass (black squares in Fig. [^, namely the deuteron binding energy and 
the singlet neutron-proton scattering length, are known accurately. Thus we abstain from a display of the sensitivity 
of those values to the uncertainty in the input. For the unphysical pion masses, the uncertainty in the input data is 
significant. For each cutoff, we thus obtain LECs not only for the central values but also for the boundaries of the 
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TABLE II. The LO LECs Cs,t and Di [GeV] for real (m^ = 140 MeV) and lattice (m^ = 510, 805 MeV) nuclei for various 
values of the momentum cutoff A [fm“^]. yields the triton as the ground (excited) state. 


rrin 

A 

Gi,o 

Go.i 

Di 

Dt 

140 

2 

-0.142 

-0.106 

0.068 

- 


4 

-0.505 

-0.435 

0.677 

- 


6 

-1.09 

-0.986 

2.65 

- 


8 

-1.90 

-1.76 

7.82 

- 

510 

2 

-0.145 

-0.130 

0.157 

-0.120 


4 

-0.438 

-0.412 

0.907 

-0.441 


6 

-0.889 

-0.853 

3.21 

-0.855 


8 

-1.50 

-1.45 

9.44 

-1.27 

805 

2 

-0.148 

-0.138 

0.071 

- 


4 

-0.405 

-0.388 

0.354 

- 


6 

-0.789 

-0.766 

1.00 

- 


8 

-1.30 

-1.27 

2.22 

- 





FIG. 2. Dependence on the cutoff A (in MeV) of the LO LECs of/PEFT: mjvA“^C'i,o (left panel), mjvA“^C'o,i (middle), and 
mivA“^Di (right). The squares (m,r=140 MeV), circles (510 MeV), and pentagons (805 MeV) represent values fitted to the 
central values of the shallowest two-nucleon S-matrix poles in the singlet and triplet channels. D\ is adjusted to the triton as 
the ground (full circles) or (for m.n-=510 MeV only) first-excited (empty circles) three-nucleon state. The shaded uncertainty 
is obtained by varying the input data within its margin of error. 


two- and three-body binding energies. In the case of I?i, we fix the two-nucleon LECs to their central values when 
varying B,^ within its error margins. The widths of the blue (to^= 805 MeV) and red/gray (510 MeV) bands in Fig.[^ 
represent how input-data uncertainty translates into coupling strength uncertainty. 

Some aspects of the cutoff dependence of the LECs shown in Fig.j^can be understood from general arguments. For 
a non-derivative four-nucleon LEC C (multiplied by A^/(167r^/^) as in Eqs. (|^ and ([^) which determines a scattering 
length o, an expansion in powers of relative momentum over A of the loop integrals appearing in the T matrix gives 

ES] 

mp,A-2c(A) =Bo + %+0 ((aA)"") , (16) 

where 9i are regulator-dependent numbers of 0(1) that depend neither on a nor on and thus also not on This 
large-A behavior is apparent in the left and middle panels of Fig. where we display TOnA^^Cs^t rather than Cs,t- 
As we can see, all curves approach a limit 0q ~ 0.7, at a rate that depends on a. The different sign of the scattering 
length in the singlet channel for physical results in a different approach to the asymptotic value compared to the 
other channels and pion masses, where relatively shallow bound states exist. 

We can also gain some insight into the cutoff dependence of Di. In the absence of a three-nucleon force, the triton 
spectrum depends sensitively on A, indicating a lack of renormalizability. The example of = 510 MeV is shown 
in Fig. When Di — 0, the open circles on the dotted line show an almost exponential dependence of the ground 
state on the cutoff. As indicated by the filled circles on another dotted line, around 1.2 GeV a second bound-state 
pole emerges, which also becomes increasingly more bound as the cutoff increases. The pattern repeats as the cutoff 
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FIG. 3. Dependence on the cutoff A (in MeV) of the three-nucleon bound-state spectrum B-j. (in MeV) in the triton channel 
at m.n-=510 MeV. The neutron-deuteron threshold is indicated by a dashed line. Empty (filled) circles mark ground (excited) 
bound states. The dotted lines are for Di = 0. For the full lines the three-nucleon force fixes the shallow state to B^,. 


increases further. Renormalizability can be achieved with the non-derivative three-body force [2m23j 

m„A-2Di(A) = F(A/A*), (17) 

where A* is an WTr-dependent parameter that determines the three-body spectrum and F is a, dimensionless function 
that depends on the regulator and on which state is kept at the observed B,^. Accordingly, in the right panel of Fig. 
we display TOnA“^£)i. 

For all values of m,r we fit Di to ensure the deep bound state remains at the observed value of in which case 
F increases monotonically with A. The resultant values of the LEG define bands which vary significantly in width 
with TO^. For m,n.=805 MeV (blue band around pentagons in the right panel of Fig. and the physical m.,^ (black 
squares), the band width is narrow relative to m,n.=510 MeV (gray band around circle^. All three bands correspond 
to repulsive interactions. This means that without a three-nucleon force there is a three-body state which is more 
bound than the observed triton. This state is then “raised” to the triton by the repulsive interaction. 

Since additional three-body bound-state poles appear at the two-body threshold with increasing cutoff, it is possible 
to renormalize Di to a shallow state instead. In this case, F changes. In the example of Fig. we can fit Di{K) so 
that the first excited state is “lowered” to the triton level as indicated by the filled circles on a full line — we label 
the corresponding values of Di as in Table [ll] and Fi g. [3} In this case, two states remain bound: the triton and a 
deep state shown by the empty circles on full line in Fig. the latter with a binding energy that goes from 170 MeV 
at A = 400 MeV to 900 MeV at A = 1.6 GeV. The increased binding of the deep state compared to its binding when 
Di = 0 shows that the force is attractive in this range of cutoffs, and it indeed has an opposite sign to the force that 
keeps the ground state at B^, as seen in Fig. 

The functional dependence of the three-nucleon LEG in Fig. is, by construction, identical to that of the binding 
energies on the cutoff when Di = 0. The latter grows faster than quadratic (upward bending of the black dotted line 
with empty circles in Fig.|^. This deviation is consistent with the increase of Di in Fig. which is not just quadratic 
but receives contributions from higher powers of A. For both fitting choices, we find the uncertainty in Di by taking 
B^ G [15.8,24.8] MeV (see Table |l]) at m,r=510 MeV. It is considerably larger when a repulsive three-nucleon force 
is used, as shown by the width of the gray band in the right panel of Fig. which is much larger than the red 
band that represents the variation in the attractive force strength. In contrast to the log-periodic behavior of the 
three-body force as a function of the cutoff found in Refs. piIEH] we find both, the central value and the uncertainty, 
to increase monotonically with A for all ttItt except for the calibration to the excited state (empty circles in Fig. [^. 
No discontinuities at critical values of A are observed because the eigenstate we chose to fit Di was always either the 
ground or the first excited state. A log-periodic F, as in Refs. I1IH13], is found if the LEG is calibrated always to the 
shallowest state in the spectrum. As a consequence, after renormalization the smallest binding energy is fixed, while 
states accrete from very large binding energies at the critical values of A. 

The significant difference in uncertainty of the three-nucleon-force parameter when fitted to the ground or excited 
state is related to the functional dependence of those states on A. In the vicinity of a critical A where an additional 
state enters the spectrum, the eigenvalue of the excited state increases much slower than that of the ground state 
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(compare slopes of the dotted lines in Fig.[^. The respective three-nucleon interaction strength inherits a larger slope 
if the ground state is fitted relative to calibrating the excited state. Since both input and regulator variation represent 
a change in the renormalization scheme, the observed difference in uncertainty is a consequence of the differences in 
slope. 

Different values of A and different regularization schemes correspond to different models of the short-distance 
behavior of the theory. These models might allow for deeply bound states in the deuteron, triton, and alpha-particle 
spectra. A tenet of EFT is that high-energy, or short-distance, phenomena can be accounted at each order by the 
most general interactions consistent with symmetries and required by RG invariance. In the specific case, we use this 
tenet to conjecture that low-energy observables, such as the nD scattering lengths, should be independent of whether 
we fit the triton to the deepest, second-deepest, ..., or shallowest state. This has been seen in simple explicit examples, 
such as Ref. [61) . where invariants of few-body spectra were analyzed with respect to changes in the short-distance 
structure of the employed models. It is not the scope of this work to assess differences between the various schemes, 
and hence we employ repulsive three-nucleon forces consistently in all calculations below. 


V. RESULTS 


The Phillips [33] and Tjon HD correlations are non-trivial features of nuclear physics. Their sensitivity with respect 
to is analyzed here. In addition, we consider the quartet three-nucleon channel which is less sensitive to the short- 
distance structure of the interaction, j.e., no three-nucleon interaction contributes up to high order. With these 
predictions, we conclude that key nuclear properties are, qualitatively, insensitive to — a conjecture based on the 
universal EFT approach. We compare the results for TO7r=5I0 MeV and m,n.=805 MeV with LO ^EFT results at the 
physical pion mass to make similarities and differences explicit. 

Besides identifying the peculiarities of large pion masses, we predict the outcomes of “experiments” in these hy¬ 
pothetical worlds. As described in the previous section, for each cutoff, i.e., model for unobservable short-distance 
structure, yPEFT differs in its coupling constants. If these models differ in predictions by a finite amount, it is this 
amount that quantifies the theoretical error. In the physical world, theoretical error estimates of this kind were used 
previously to make predictions through correlations (see e.g. Refs. [151135]). Since the theoretical error for and 
Rq, is large relative to the experimental one, the Phillips and Tjon lines at LO in /EFT do not constrain observables 
further at physical m^n-. At larger however, the lattice uncertainty is still significant (see Table |T]) and, a priori, 
there is no reason why yPEFT should not be able to constrain observables more tightly through those correlations 
than solely by the “experimental” error. 

We assess sensitivity of results to higher-order terms in the EFT expansion by a variation of the cutoff-regulator 
parameter in the range A G {2,4,6,8} fm“^. This range includes the critical value for appearance at the physical pion 
mass of an excited state, when Di — 0 (see discussion of Fig. in Sec. IV). For lattice pion masses, the upper limit 
is 2-3 times the expected breakdown scale, where in general we see signs of convergence in observables. Although we 
might ideally want even higher cutoffs at the expense of further computational time, our estimate of the truncation 
error is probably not an underestimate because we include cutoff values below the expected breakdown scale. Such 
low cutoffs introduce larger errors than the truncation. A more reliable error estimate has to await higher-order 
calculations where the breakdown scale manifests itself. 


A. The Three-Body Sector 

The physical nucleon-deuteron system splits into two significantly different spin channels: doublet (or triton) with 
s = 1/2 and quartet with s = 3/2. The former (latter) supports (does not support) a bound state. In the doublet 
channel, an additional counterterm enters at LO — Di term in Eq. (§ — while the quartet channel is renormalized 
with Cs,T, only. The consequences to large TOtt are the subject of the following two sections. Since we include no 
Coulomb interactions, our results at physical pion mass apply only to neutron-deuteron scattering. 


1. Neutron-deuteron ^ 83/2 channel 

The phase shifts in the quartet channel can be calculated in LO solely on the basis of two-nucleon input. In Fig. [^ 
we show our phase shift results for elastic nD scattering below 10 MeV, calculated with the MGM. For all three pion 
masses, the cutoff variation between 2 fm“^ and 8 fm“^ is shown by green (m7r=140 MeV), red (TO.n.=510 MeV), and 
blue (to^= 805 MeV) bands. The upper (lower) edge corresponds to 8 (2) fm“^. For the physical m,rj we compare 
our results to previous LO and N^LO /EFT calculations [TT] [22] (black solid and dashed-dotted lines) obtained from 
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FIG. 4. Elastic neutron-deuteron scattering phase shift 5 in the spin-quartet * 83/2 channel (in degrees) for various pion 
masses as function of the center-of-mass energy Ec.m. (in MeV). The green (m.n- = 140 MeV), red (m^=510 MeV), and blue 
(m^=805 MeV) shaded areas are the LO/PEFT results of the RGM for a cutoff in the range [2,8] fm“^. For the overlapping 
m,r = 510, 805 MeV uncertainty bands, the upper (lower) edge, corresponding to A = 8 (2) fm“^, is given by dashed lines. For 
m,r=140 MeV, the solid (dashed-dotted) black line represents LO (N^LO) /fEFT results from the solution of the STM equation 

[niiii]- 



FIG. 5. Dependence on the cutoff A (in MeV) of the quartet C^OnD, left) and doublet (^Uno, center) neutron-deuteron 
scattering lengths (in fm), and of the alpha-particle binding energy (right, in MeV), for rtir: = 140 MeV (black squares), 
510 MeV (red circles), and 805 MeV (blue pentagons). 


the solution of the Skorniakov-Ter-Martirosian (STM) equations. The difference between these curves is, of course, 
a good reflection of the uncertainty of the LO calculation at the physical pion mass. Our result has the correct 
energy dependence and lies between the two curves. Our error band accounts for cutoff variation but not numerical 
uncertainty. The latter is included in the postdiction of the nD scattering length given below. The energy dependence 
and band widths are similar for the three values of the pion mass we consider. This suggests an invariance with 
respect to m,n. of the uncertainty — and therefore the convergence rate of the EFT. 

We extract a scattering length at „ = 0.001 MeV. The cutoff dependence is illustrated in the left panel of Fig. 
For all values of the pion mass we observe a nice convergence pattern. Our final values are shown in the first row of 


Table III The errors are the sum of the sensitivity to higher-order effects assessed with the cutoff variation, and the 
numerical uncertainty, which we measured to be less than 1 fm. They are of similar size for the three pion masses, as 
for the phase shifts at higher energies. 

The quartet scattering length is an example of what is sometimes called a low-energy theorem: to a high order 
it is entirely determined by LECs fixed in other processes. The value we obtain for at physical pion mass is 
consistent with the -^EFT postdictions [20] of "‘onD = 5.1 ± 0.80 at LO and '^OnD = 6.4 ± 0.020 at N^LO, and with the 
experimental value [79]. We find a slow decrease with rriTr, but no significant change, which could have arisen if there 
were a shallow bound state in this channel. The ERE should apply below the deuteron break-up threshold, where 
the deuteron can be treated as a single body. We might expect that, barring some fine-tuning, the size of the ERE 


parameters is set by the deuteron break-up threshold, fcpn — The numbers in Table III show indeed 
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TABLE III. Leading-order postdictions (ytEFT) and predictions (?fEFT) for the quartet and doublet neutron-deuteron scattering 
lengths '^flnD and ^OnD at three pion masses, in comparison with experiment and LQCD. The theoretical uncertainty considers 
cutoff variation between 2 fm“^ and 8 fm“^, model-space truncation, and LQCD-input variation. 



tIEFT 

^EFT 

[MeV] 

140 

510 

805 

"‘onD [fm] 

5.5 ± 1.3 

2.3 ± 1.3 

1.6± 1.3 

[fm] 

0.61 ±0.50 

2.2 ±2.1 

0.62 ± 1.0 


experiment [79] 

LQCD 

'‘onD [fm] 

6.4 ±0.020 

7 

7 

[fm] 

0.65 ± 0.040 

7 

7 


very good agreement with the expectation I^Onol = C>(l/fcpn)- 

Both the convergence pattern with the cutoff (reflected in error bands) and the natural size of the resulting ob¬ 
servables suggest that yt'EFT behaves in similar ways to ^EFT. There is no evidence that observables in this channel 
require a different treatment from an EFT point of view for the larger pion masses, i.e., the same power counting is 
applicable. 


£. Neutron-deuteron ^Si /2 channel 


As precise calculations of the three-nucleon system became possible in the late 1960s, correlations were observed 
among certain three-body observables calculated with a variety of potentials fitted to two-nucleon data. The best- 
known example is the Phillips correlation [33] between the doublet scattering length and the triton binding 

energy B.^.. In ^EFT these correlations are understood by the fact that, if the two-body input is fixed, three- 

body predictions in the doublet channel still depend on one parameter in LO, which determines the three-body force 
in Eq. 0. As this one parameter is varied, three-body observables sensitive to the LO three-body force all change 
in a correlated way. t^EFT, by construction, predicts the correlations, and we establish here similar higher-order 
uncertainties for all pion masses. 

Fixing the LO three-body parameter to one datum, other observables are calculated as for the quartet channel. 
The cutoff dependence of the doublet neutron-deuteron scattering length is shown in the middle panel of Fig. [^in 
the case where the three-body parameter is determined by B^, as described in Sec. IV Again signs of convergence 
are visible, but not as clearly as for the quartet scattering length. The lowest cutoff of 2 fm“^ is not clearly above 
the breakdown scale at m,r = 805 MeV, and indeed it generates significant errors. Therefore, for this pion mass we 
consider only cutoffs 4 fm or higher in the following analysis. 

One also expects the values of the doublet scattering length to be correlated with the triton binding energy. This 
is particularly clear when \B.j, — Bj^\ B^, in which case the triton can be described as a neutron-deuteron bound 
state and the small binding translates into l/^OnO ^ 1- But this correlation is observed also beyond this region: in 
Fig.i the Phillips correlation is shown for the three pion masses. As the renormalization condition fixes the binding 
energy but does not eliminate a residual cutoff dependence, which can only be removed by higher-order interactions, 
the correlation is manifest as a band of finite width rather than a one-dimensional line. This width represents the 
theoretical error at LOX'EFT. The bands were mapped out by a line for each cutoff A € {2,4,6, 8} fm“^. Each such 
line is parametrized by a factor multiplying the three-body interaction. At TO7r=805 MeV, the correlation is about to 
break down for the lowest cutoff, A = 2 fm“^ (blue dashed line in Fig. [^, which is another evidence that this value 
cannot be considered representative of the EFT truncation error. 

We extract the values for ^OnD shown in the second row of Table III 


In the doublet channel, a too-small model 
space can be more easily detected than in the quartet channel from an under- or over-bound triton. As a consequence, 
the numerical KGM uncertainty is about 0.1 fm and therefore small relative to higher-order effects which are taken 
as the width of the band: 0.26 fm for TO,n.=510 MeV, and 0.13 fm for TO7r=805 MeV. The approximately constant 
width of the Phillips band for all three suggests invariant higher-order uncertainty with increasing Br,,. Since 
higher-order effects scale with momenta — those of nucleons increase in the triton as B^ increases — the band should 
intuitively widen towards larger B^. In effect, lattice input uncertainty dominates the theoretical error. For the two 
unphysical pion masses (center and right panels), the gray-shaded areas represent data uncertainty given in Table [I] 
The intersections of the edges of the error bands with the correlation bands define areas (gray areas bounded by 
horizontal dashed lines) in the B^ —^ o„d plane which contains pairs of values that are consistent with all other data 


points. The total theoretical uncertainty as estimated in Table III includes the error in the LQCD input, 1.71 fm for 
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FIG. 6. Correlation between the doublet neutron-deuteron scattering length (in fm) and the three-nucleon binding energy 
(in MeV). The green (m.,r=140 MeV, left panel) and red (m.,r=510 MeV, center) shaded areas are the LO^EFT results of 
the RGM for a cutoff in the interval [2, 8] fm“^. The blue (m.,r=805 MeV, right panel) shaded area is the LO^EFT results of 
the RGM for a cutoff in the interval [4,8] fni~^, with the cutoff A = 2 fm“^ shown as a blue dashed line. For m,r = 140 MeV 
(left panel), experimental data are marked with a red dot and blue dashed lines display LO tIEFT results obtained with the 
STM equation for sharp cutoffs of 140 and 900 MeV |16| . The gray shaded area (center and right panels) marks the lattice 
uncertainty in Values for ^Ono between the horizontal dashed lines are consistent with all other low-energy data. 


niTr = 510 MeV and 0.26 fm for m,n. = 805 MeV. 

The calculation at physical ttItt (left panel in Fig.j^ serves as a benchmark. The dashed lines represent the solution 
of the STM equation at LO in /EFT when a sharp cutoff was varied between from 140 to 900 MeV [T^ . Experiment is 
represented by the red dot. Our band is consistent with both. The small scattering length compared to the break-up 
threshold inverse momentum, 1/fcpn — 2.2 fm, is a sign of a zero of the T matrix near the nD threshold. As discussed 
in Ref. [55], such a zero is located at ~ ^pn with respect to the origin of the complex relative-momentum 

plane, when [fcol ^ ^pn and the deuteron can be treated as elementary. A consequence is a large effective range 
^?"nD ~ —(^OnD^o)”^ a Small radius of convergence of the usual ERE. Data suggests |fco| ^ 20 MeV on the 
imaginary momentum axis, and indeed this is what was found by explicit calculation already many years ago [80j . 
The negative slope of the Phillips line indicates that, as the three-body force is changed so that the triton gets more 
bound, this pole crosses threshold. The zero remains in the region of validity of the elementary-deuteron theory for 
1 or 2 MeV around the physical value of . In this region a modified ERE [68] holds [66] ■ 

In the center and right panels of Fig. we extend LQCD to the realm of few-body scattering, which is not as easily 
accessible directly on the lattice. The negative slope of the Phillips line persists but the line moves up in the ^ OnD 
plane, and it gets flatter, as TOtt increases. The “measured” triton binding energy rises monotonously (Table |l]), the 
doublet scattering length first increases then decreases (Table [TTl| ) . For = 510, 805 MeV, the accidental zero of the 
nD scattering amplitude is no longer clearly present in the region where the deuteron can be taken as elementary, 
and no particularly large effective range is expected. 

The increasing with increasing B^ from 140 MeV to 510 MeV pion mass is opposite to the trend found for 
fixed pion mass, as identified above. It is not the triton binding energy but the triton-to-deuteron binding ratio 
shown in Table |l] which decreases with increasing ^OnD- This ratio is important because it measures the motion of 
the “experimental” point in the B^ —^ UnD plane: B^ influences (together with B^^) the position of the line, and B,^ 
fixes the position along the line. As the ratio B^/B^ decreases from m,r=140 MeV to m,r=510 MeV and increases 
from 510 MeV to 805 MeV pion mass, the three-nucleon bound state comes closer to and farther away from the nD 
threshold. 

In particular, the relatively large scattering length at m,r = 510 MeV reflects a less-bound triton relative to the 
neutron-deuteron threshold. Once the errors are considered, B^ is just a few of MeV away from B^ at = 510 
MeV. In fact, in contrast to = 140, 805 MeV, the binding momentum of the last nucleon can be smaller than the 
deuteron break-up momentum, Khd — — B^)/3 < fcpn, in which case the ERE is likely to apply. This leads 

to the prediction 


^nD — 


^nD 


^nD ^nD 


(18) 


The hrst term gives, for the central values of the binding energies, ^OnD = 1-6 fm with a correction of about 50% from 
the second term if prnDl l/^pn- This estimate agrees with the central value calculated with the full three-body 
dynamics given in Table [ITT| Barring signiheant shifts in the central values as lattice errors are reduced, in this lattice 
world the triton can be viewed as a two-body halo system. 
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FIG. 7. Correlation between the three- (43 t, in MeV) and four- {B^, in MeV) nucleon binding energies. The green 
(m,r=140 MeV, left panel), red (m.,r=510 MeV, center panel) and blue (m.,r=805 MeV, right) shaded areas are the RGM- 
LO^EFT results for a cutoff in the interval [2,8] fm“'^. For m^=140 MeV (left panel), the green uncertainty band represents 
sensitivity to the cutoff and to the renormalization input (whether or ^Unp). Experimental data are marked with a red 
dot, and the blue dotted (dashed) line represents LO ylEFT results from Ref. [2^ using ^anp,^anp (^anp,BD) as input. The 
gray shaded areas in the center and right panels mark lattice uncertainty in and . Values for B^ between the horizontal 
dashed lines are consistent with all other low-energy data. 


It is an open question whether there is a pion mass, possibly around 510 MeV, where the triton converges to the 
deuteron threshold. If there is, we would be witnessing a qualitatively new phenomenon in few-nucleon physics. The 
Efimov effect in the three-body system is a prominent example of universal feature emergent from the unitary limit 
in the two-body sector. A pion mass which produces the analog three-body unitarity, l/^OnD —> 0, would be a world 
where the four-body system exhibits an Efimov-type spectrum. 


B. The Four-Body Sector 

While there is no lattice data on three-nucleon scattering observables and thus the results presented in the previous 
subsection remain to be verified “experimentally”, i.e., with a direct LQCD calculation, there is data on the ground- 
state energy of the four-nucleon system. In this section, we find the three- and four-nucleon ground-state energies 
correlated for all three m,r- At the physical the relation is known as the Tjon line El which can again be 
explained by a variation in the single LO three-body force parameter. 

In Fig. the correlations between the ground-state energies of the three- and four-nucleon systems are shown. 
The different graphs represent results for the three pion masses. We observe an increase in alpha-particle binding 
in step with the increase in triton binding energy, which is not surprising in pionless EFT because with fixed two- 
nucleon input it is the same three-body force that controls the binding of the three- and four-nucleon systems. The 
correlation is manifest in a band, not a line, and the width of the band measures the theoretical uncertainty assessed 
via cutoff variation. With the central value of as input in the three-body force, the dependence on the cutoff of 
the alpha-particle binding energy is shown on the right panel of Fig. The slope of the correlation lines — as 
before, each line is parametrized by a variation of the Di three-body interaction strength — further affects the LO 
ydEFT uncertainty. The larger the slope, the larger the uncertainty in B^ due to the uncertainty in B^. 

For the physical pion mass, our error band does not agree well with the LO results of Ref. [29] shown in the left panel 
of Fig.[^ In Ref. |5S] the alpha-particle binding energy was found by a solution of the Faddeev-Yakubovski integral 
equations with a Gaussian regulator on the relative incoming (p') and outgoing (p) momenta, exp[—(p^-|-p'^)/A^]. The 
uncertainty was assessed by a cutoff variation A G [8,10] fm“^, thus excluding a reported stronger cutoff dependence 
for A < 8 fm~^. The cutoff variation was then deemed small compared to the higher-order uncertainty estimated 
by changing the two-body input: the two curves obtained from B^ and are represented by the blue lines in 
Fig. 0 We have similarly examined the input dependence: for A = 8 fm we replaced with ^Onp and found B^ 
(upper bound of the correlation band for to^= 140 MeV, green area, left panel in Fig. larger by 2 MeV. Even with 
this extended variation of the renormalization scheme, the two uncertainty bands do not overlap. In contrast, the 
current KGM results for the Tjon correlation band are consistent with the previous KGM LO-^EFT calculation of 
Ref. [35]. The convergence of B^ to the physical value when the NLO potential is iterated [35] suggests that in both 
LO calculations the theoretical error as shown by the band widths in Fig. [^is a lower bound. For our theoretical 
error estimates, we interpret RGM and Faddeev calculations, i.e., different regularizations and model-space cutoffs, 
as different renormalization schemes. For B^ and physical m^^., the uncertainty is thus given by the spread of results 
of both methods (difference between short-dashed blue line and lower edge of the green band). 














16 


TABLE IV. Predictions for the four-nucleon binding energy and the universal alpha-to-triton ratio B^jB^ from LO 
pionless EFT at three pion masses, in comparison with experiment and LQCD. The theoretical errors include numerical and 
EFT uncertainty. The uncertainty in the fractions (lines 2 and 4) adds independent errors in quadrature. 


niTT [MeV] 

140 

510 

805 


/EFT 

/PEFT 

B^ [MeV] 

24.9 ±4.3 

35 ±22 

94 ±45 


2.9 ±0.51 

1.7±1.1 1.8±0.9 


experiment 

LQCD 

B^ [MeV] 

28.3 

43.0 ± 14.4 107.0 ±24.2 


3.34 

2.1 ±0.85 

2.0 ±0.6 


For unphysical pion masses we indicate, as before, the uncertainty in LQCD energies by gray bands in the center 
and right panels of Fig. a vertical band for and a horizontal band for B^. Values for B^ in an interval 
bounded by the intersection of the upper (lower) edge of the ydEFT uncertainty band with the right (left) boundary 
of the band of LQCD-allowed B^ values are indicated by horizontal dashed lines. This range is slightly larger than 
the constraint already given by “experiment” for 771,^=805 MeV, and slightly narrower for m,n.=510 MeV. However, 
given the renormalization-input dependence seen at physical pion mass, we estimate the theoretical uncertainty by 
conservatively doubling the width of the KGM correlation band, plus 2 MeV as an upper bound of the numerical 
uncertainty (see Fig. for this estimate), plus the experimental LQCD uncertainties in B^. 

The results for the alpha-particle binding energy are summarized in Table |IV| The predicted value is taken as the 
central value in the uncertainty band. The EFT results—absolute binding energies and the ratios—are consistent 
with experiment at physical pion mass and with LQCD at higher masses given the uncertainty estimates on both, 
experimental and theoretical side. One should keep in mind that the experimental number reflects the additional 
effects of the Coulomb interaction between protons, which does not enter the LQCD results. As discussed in Sect. |TTj 
Coulomb effects should be of higher order in the relatively tight helion and alpha-particle ground states. 

As expected, the B^ — B^ correlation has a positive slope for any pion mass. For each correlation band, we define 
the slope with a linear regression through the B ^, B^ pairs predicted withyf'EFT for all cutoff values which are within 
data uncertainty (gray areas) only. At physical our calculation yields a smaller slope (« 3.6) than Ref. [29] 
(« 3.8). With increasing m,r, the slope is found to decrease, « 2.1 for m,r=510 MeV and m7r=805 MeV. In other 
words, the ratio B^jB^ does change with consistent with the lattice extractions, as shown in Table IV 


Since in obtaining the Tjon line the three-body force is being varied, the structure of the line (slope, curvature, 
intercept) must depend on the two-nucleon interactions. Indeed, from the ratios listed in Table |IV| we infer that 
whatever leads to the different ratios between the deeper two-nucleon state — recall that for unphysical masses, the 
interaction in the channel sustains a bound state, see Table |T]— and the triton is not the main factor behind the 
change in the slope of the Tjon line. As for the Phillips line, the structure of the Tjon line depends on both pieces of 
two-nucleon input. For example, in Fig. [^ we demonstrate the sensitivity of the slope of the Tjon line with respect 
to the pole in the spin-singlet two-nucleon amplitude. The triplet binding energy, i.e., the deuteron was held fixed to 
the lattice central value at m7r=510 MeV , B^ = 11.5 MeV. Three cases are shown for A = 4 fm“^, corresponding 
to different calibrations of Cop, the LEC controlling the channel: i) a singlet neutron-proton state with binding 
energy of approximately 11.5 MeV, i.e., the deuteron energy; ii) a shallow bound singlet state of B^ « 0.5 MeV; 
and Hi) an unbound singlet state. Within the considered range for B.^. a linear regression to the dependence of B^ 
on B.J, is appropriate. When the singlet two-nucleon state is very close to threshold, the slope is found maximal, 
AB^/AB.J, « 3.0 (red dashed line in Fig. |^. If the interaction is tuned away from this critical point, either to 
produce no bound singlet (red dotted line), or a state with identical binding energy to the triplet (red solid line), the 
slope parameter decreases. Naively, one might have expected a monotonic dependence of the slope on the strength 
of the two-body attraction. A larger two-body attraction requires a more repulsive three-body force to fix the triton. 
The contribution of this extra repulsion should be stronger in the four-particle system and hence the latter should 
not be as deeply bound. The non-monotonicity found above remains to be explained in a more general context taking 
into account the possibility of a four-body Efimov effect mentioned at the end of Sec. V A 2[ 


VI. CONCLUSION 

We have adapted pionless effective field theory, /EFT, to describe LQCD data at unphysical pion masses, dubbing 
it yf'EFT. For the first time predictions were made for a nuclear reaction, nucleon-deuteron scattering, in lattice 
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FIG. 8. Correlation between the three- (Bt, in MeV) and four- in MeV) nucleon binding energies (Tjon line) for three 
values of the binding energy of the singlet two-nucleon state: same as the deuteron binding energy (continuous line, slope 
~ 1.9); shallower than the deuteron (dashed line, slope « 3.0); and unbound (dashed line, slope ~ 2.4). Results were obtained 
for ten values of the three-body interaction strength at A = 4 fm“^ for = 11-5 MeV (m^=510 MeV). 


worlds where the pion mass is 510 and 805 MeV. Furthermore, the Phillips and Tjon correlations were obtained at 
these high pion masses with leading-order uncertainties of similar size and thus offer no indication of a significantly 
different convergence rate of the respective EFT expansions. The alpha-particle binding energy was found in good 
agreement with direct lattice measurements, which reassures us of the applicability of ydEFT. It also strengthens 
confidence in the LQCD results [37l [38] themselves, despite the apparent subtlety in identifying energy plateaus in 
the data. 

Our work thus suggests that yP'EFT applies to light nuclei independently of the exact LQCD data used as input. 
The calculations presented here could be repeated if those values change or if new values of the pion mass are explored. 
While this manuscript was being written, new data have appeared for TOtt = 300 MeV m, which do not quite fit with 
the trend of increasing binding with pion mass but show a pattern of dependence on A similar to the one found at 
higher pion masses [sangQ More problematic is that another lattice collaboration does not find bound states 
over a wide range of pion masses. Because the latter lattice data are processed through an (unobservable) potential, 
uncontrolled errors are introduced. Still, it is prudent to consider the specihc numbers available from LQCD so far 
as illustrative only, and focus instead on the qualitative insights they bring into nuclear physics [2|. 

While much of the underlying structure of light nuclei seems to remain the same at unphysical pion mass, existing 
LQCD data give some hints of subtle changes. In the studied lattice worlds, the triton and alpha-particle binding 
energies are larger than in the real world, but their ratios to the deuteron binding energy become smaller. In contrast 
to the quartet neutron-deuteron channel, where we detected no qualitative changes, the accidental zero of the doublet 
T matrix that exists for physical pion mass near threshold seems to disappear. It is replaced by effective-range 
parameters that suggest a more prominent neutron-deuteron halo character for the triton. LQCD data with smaller 
errors, and at other values of the pion mass, would allow firmer conclusions about the organization of nucleons in the 
triton and its implications for the alpha particle. 

In an upcoming project, NLO calculations will assess the convergence rate and breakdown scale of the EFT. They 
will also reduce some of the uncertainties in the extrapolation of LQCD data. In the longer term, the application of 
X'EFT to systems with more than four nucleons could guide the lattice effort to the relevant few-body observables to 
be measured in order to pin down additional LECs needed to understand nuclear structure. 


^ Note that the central value s fo r deuteron and triton energies m suggest a triton with a last nucleon that is even less bound than at 
rrin = 510 MeV. From Eq. |l8| we then expect — 2 fm with a correction from the effective range of perhaps 40%. Alas, in this 

lattice world, too, the large lattice errors do not yet allow firm conclusions about the two-body halo character of triton. 
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